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Abstract: The effective four-dimensional, linearised gravity for a brane world model 
with higher order curvature terms and a bulk scalar field is analysed. Large and small 
distance gravitational laws are derived. The model has a single brane embedded in a five- 
dimensional bulk spacetime, and the scalar field represents the dilaton or a moduli field. 
The quadratic, Gauss-Bonnet curvature term (and corresponding higher kinetic terms for 
the scalar) is also included in the bulk action. It is particularly natural to include such 
terms in a brane world model. Boundary terms and junction conditions for the higher 
order terms are given. 

The extra terms allow additional solutions of the field equations, which give better agree- 
ment with observational constraints. Brans-Dicke gravity is obtained on the brane. The 
scalar and tensor perturbations are affected differently by the higher gravity terms, and this 
provides a way for the scalar modes to be suppressed relative to the tensor ones. Another 
new (but less useful) feature is the appearance of instabilities for some parameter ranges. 
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1. Introduction 

There has been considerable interest in the idea Q that we may live on 3+1 dimensional 
brane embedded in a 4+1 dimensional bulk spacetime. Despite the extra dimension, it 
is still possible for the effective gravitational theory on the brane to closely resemble that 
which is observed in our universe. The models we will be considering in this paper resemble 
the well known Randall-Sundrum II brane world model [Q]. In this scenario the graviton 
spectrum has a zero mode which, due to the warping of the bulk space, is localised on the 
brane. This is interpreted as the four-dimensional graviton. At large distance scales the 
effects of the other modes are suppressed, and something close to four-dimensional general 
relativity is obtained for brane-based observers Gravity will have unconventional be- 
haviour at short distances, and so to be compatible with gravity experiments, 'short' needs 
to of order 0.1mm or less. 

The brane world scenario is loosely motivated by string theory. However, if such a 
model were actually derived from string theory, the bulk space would contain not only 
gravity, but many other fields as well. The simplest example of this would be a scalar 
field such as the dilaton, or a moduli field coming from the compactification of other extra 
dimensions. One of the purposes of this paper is to investigate the implications of bulk 
scalar fields for brane world gravity. For simplicity I will include just one such field, and 
consider only a special class of solutions for which the field equations can be easily solved. 

Scalar-tensor gravity has previously been studied in a brane world context [Q], although 
most previous work has considered only the graviton modes. However, given the strong 
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constraints placed on scalar modes by solar system measurements , we should not ignore 
them if we are studying a brane world which is supposed to represent our universe. 

Another natural extension of the usual brane world model is to add higher order cur- 
vature terms to the action. One such term is the quadratic Gauss-Bonnet (or Lanczos Q) 
term 

£gb = R^- ^RabR"' + R'''""'Rabcd . (1.1) 

Like the Einstein-Hilbert action this gives rise to field equations which are divergence free, 
and do not contain any third or higher order derivatives of the metric. They therefore 
have all the desirable properties of the Einstein equation (such as a ghost-free vacuum and 
conservation of energy momentum) . Note that this is the only quadratic curvature term 
which has all these properties. In four dimensions its contribution to the field equations is 
trivial, and so it is usually ignored. 

The implications of this term for brane world gravity have been discussed in the lit- 
erature js], |l^, |ll|, 12 1 . Note that the derivation of the brane junction conditions is 



more complicated than in the conventional brane world ^ , 13 1 , and this led to some subtle 
mistakes in early works 

Like the branes themselves, higher order curvature terms (and corresponding higher 
order scalar kinetic terms) are also motivated by string theory. It is therefore very natural 
to consider them in brane world models. For simplicity, we will only consider the Gauss- 



Bonnet term and the corresponding scalar field terms in this paper (see ref. |1C] for a 
discussion of even higher order curvature terms in a brane world context). 

The layout of the paper is as follows: In section ^ the details of the model are discussed. 
Solutions of the field equations are described, with and without the higher order terms in the 
action. Not surprisingly, the addition of the higher order terms produces extra solutions 
to the field equations. For want of a better name we will call these the 'Gauss-Bonnet 
branches', and the other solutions the 'Einstein branch'. Linearised brane world gravity 
is reviewed in section |^. Linear perturbations of the 'Einstein' solutions are considered 
in section and the more interesting and better behaved 'Gauss-Bonnet' solutions are 
analysed in section |5|. Several of the solutions are shown to be unstable, although some 
of these instabilities can be removed by adding appropriate terms to the brane part of the 
action. These stabilised solutions are discussed in section |^. The last section contains a 
summary of the results. The bulk field equations and junction conditions (together with 
corresponding boundary terms) are listed in the appendix. 

2. Bulk Solutions 

We will consider a five-dimensional toy model with a single, Z2-symmetric brane and one 
bulk scalar field, <I>, which is conformally coupled to gravity. We will be working in the 
Jordan frame (in a string theory context, and when <I> is the dilaton, this is known as the 
string frame). We start by considering the first order action 

J d'x^e-^" {R - 4cu(V$)2 - 2A} - J d^x^/^e'^* {2K + T} , 

(2.1) 



2 
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where hab = gab — nafib is the induced metric on the brane and Kab = h'^o^c'nb is the 
extrinsic curvature. M is the five-dimensional Planck mass, A is the bulk cosmological 
constant, and T is the brane tension. The brane is being treated as a boundary of the two 
halves of the bulk spacetime, and the Gibbons-Hawking boundary term has been included 
to give a consistent action |14]. In this paper we take the brane normal to be pointing 



away from the brane on both sides of the brane. This is an unconventional choice, but it 
is easier to visualise. 

The general second order action is 

S2 = Y J d^x^e-^'^ {ci£gb - 16c2G«6V"W''$ + IGc^iV^ fv"^^ - 16c4(V$)^} 

-M j /:E^/^e-2*|^ci {sKK^bK^' - 2K\K''Kac - - GG'^'Kab) 

16 1 

- 16C2(,K,b - Khab)D''^D'^ + yC3 {{dn<^>f + SOnHD^f) | , (2.2) 

where 9„ = n"'da- The caret denotes tensors with correspond to the induced metric hab, 
and Da is the corresponding covariant derivative. 

When the scalar field is included there are a total of four possible second order terms 
which produce suitable field equations (no higher than second derivatives, etc.). For the 
action to be consistent, we must include the Gauss-Bonnet equivalent of the Gibbons- 



Hawking boundary term |13, 15], as well as the corresponding $ dependent terms (see 
appendix) . 

We see that the boundary parts of the above actions include terms which are first and 
third order in derivatives. It is therefore natural to also include an induced gravity action 

M2 



y"d^x\/^e-2* |6i^ + 62(^^)^} , (2.3) 
which is second order in derivatives. Although we will consider this term, we will be mainly 
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interested the implications of the other parts of the action, (2.1) and (p.2|). 

Scalar-tensor theories of gravity are also studied in the Einstein frame, which is related 
to the Jordan frame by the conformal transformation. The transformation is chosen so 
that ^-dependent factor in front of Einstein-Hilbert part of the action is cancelled. For 
a brane world this is either gab — > e^^^^Oab or gab e.'^^dab: depending on whether we 
are considering the bulk ( p.l[ ) or brane ( p.3| ) part of the action. So we see that there are 
actually two Einstein frames for a brane world. 

The values of the coefficients to, Ci and hi will be determined by the fundamental 
theory from which the actions ( p.l| - |2.3| ) are derived. One of the simplest possibilities is 
to suppose that $ is a moduli field arising from a toroidal compactification of a higher 
dimensional space. Consider a (5 -|- A^)-dimensional theory, whose bulk action contains 
only the Einstein-Hilbert and Gauss-Bonnet terms. Take the metric to have the form 

dsl+N = gabix)dx''dx' + e-^'^^^y^VABdX^dX'' (2.4) 

and then compactify the last dimensions (we have made the simplifying assumption 
that these extra dimensions only have one degree of freedom) [l^. Comparison of field 
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equations reveals that such a theory wih be equivalent to that obtained from the above 
actions if a; = — 1 + A''"^, 

C2 = —uJCi , C3 = {2uJ + l)ujci , C4 = —(2a; + l)a;^ci (2-5) 

and 62 = — 4cl'5i. In the — > 00 limit, $ can be interpreted as the dilaton for a theory with 
a high degree of symmetry [^]. For the rest of this paper we will assume the above relation 
between q, bi and u, and take > w > — 1. For later convenience we define a = ci/M^ 
and P = bi/M. Both of these parameters will be taken to be positive. 

Combining all the above contributions to the action ( |2.1| - |2.3| ), and including some 
brane matter gives the total action 

S = Si + S2 + 5'ind - j d'^xV^Craa.t ■ (2-6) 

The brane matter Lagrangian £mat will be treated as a small perturbation to the back- 
ground solution. For now we will set it to zero. 

In this paper we will be interested in perturbations of non-singular brane world solu- 
tions of the form 

^ g-^'^l^ldx^dx^r/^^ + dz^ , (2.7) 

with ^ = (j)Q — (t\z\. The single brane is located at z = 0, and the bulk is taken to be Z2- 
symmetric. Solutions of this form are the simplest extension of the usual Randall-Sundrum 
model Previous work []l6| , 0, |l^ on such solutions in Gauss-Bonnet brane worlds has 
used different coefficients (q) for the terms in bulk action (|2.2|), although this did not lead 



to qualitatively different results. Other, more general classes of solutions exist |12, 17], but 
we will not consider them here. 

There are problems with the above solution if o" < 0. Although (in the Jordan frame) 
there are no curvature singularities, the coupling of matter to gravity (which is proportional 
to e*) will be divergent as 2; — > 00. If A: > the time taken (as observed on the brane) for a 
massless particle to travel from this singularity to the brane will be infinite. However it is 
still possible for signals from points arbitrarily close to the singularity (and with arbitrarily 
high coupling) to reach the brane in finite time. The low energy action we are using will not 
be valid near the singularity, and so the physics of our brane universe would be dominated 
by effects that we have not included in our model. In the Einstein frames (bulk and brane), 
the curvature at this point becomes singular, and furthermore the singularity will then be 
at a finite proper distance from the brane. Therefore, in order to have a consistent single 
brane model, we will only consider solutions with o" > 0. e* is then bounded above by e*^", 
its value on the brane. 



Considering only the lowest order action (2J.), the bulk field equations (see appendix) 
are solved by 

k = -2{uj + l)a. (2.8) 

The value of a will be related to the bulk cosmological constant. We see that (for the 
a > 0, uj > —1 solutions we are considering) the sign of warp factor is the opposite of that 
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in the usual brane world model. As we will see in section ^ the negative warp factor does 
not give a suitable effective theory of gravity on the brane. 

When the higher order terms ( |2.2| ) are included, the above solution is still valid, and 
(if w / 0) two new branches of solutions appear 



(if w = then k"^ = l/[12a], the field equations are degenerate and ^ is undetermined). 
It is possible for these solutions to have positive warp factor, suggesting that a viable 
effective gravitational theory could be obtained on the brane. These new solutions will 
be referred to as the 'Gauss-Bonnet branches', and the other solution as the 'Einstein 
branch'. It is perhaps worth noting that for a more general choice of the coefficients in 
the action (|2.2|), the 'Einstein' solution will not still be valid when the higher order gravity 



terms are switched on (see e.g. [^, 18]). 
3. Linearised Brane Gravity 

We will start by reviewing linearised gravity for a brane world without scalar fields. By 
taking i2mat to be small, a perturbative analysis can be used to determine the approximate 
effective gravitational law on the brane. We will use a gauge in which the brane is kept 
straight, even when matter is present [19, [20|]. 



The general perturbed metric can be written as 

ds^ = e~'^^\'\r]^y + -f^^)dx^dx'' + 2Vf,dx^dz + (1 + ^p)dz'^ . (3.1) 

The quantities 7^jy, and ijj are all small. The perturbed brane normal is then = 
5^ sign(2:) (1 + tp/2). We use the metric convention rj^i, = diag(— 1, 1, 1, 1). 

The lapse function, ip, and the shift vector, v^, are determined by solving the com- 
ponents of the field equations which are perpendicular to the brane. If the brane does 
not bend, as we are requiring, these will not depend on the brane matter. For the Gauss- 
Bonnet-Randall-Sundrum model {k > 0) they are solved (for z > 0) by 

V' = -4^5,7 (3.2) 

Vf, = -^d^j + Bf, (3.3) 
d^%u = (3.4) 

where = Ifiu - (l/4)7r?^i/, 7 = 'n'^''7^lu, and satisfies S'^S^ = and = 0. We 

can use gauge freedom to set = 0. 

For the above choice of gauge the remaining bulk field equations reduce to 

(1 - 4aA;2) (^Sf - 4kd, + e^'^'B^^ %^ = (3.5) 

with ^4 = -q^^d^di,. There is no bulk equation for 7. The remaining gauge freedom can 
used to choose its behaviour there (but not its value on the brane) [l9|] . 
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Note that the above equation is free from third or higher order derivatives of 7^1/, and 
is similar to hnearised Einstein gravity in this respect. This follows from the fact that the 
Gauss-Bonnet term was specifically constructed to give a Einstein-like theory of gravity. 
If 1 — Aak'^ < the j^^, kinetic term in the linearised effective action will be negative, and 
so the theory will have bulk ghosts |1C, 12|. Therefore if the theory is to be stable on the 



quantum level, k must be less than l/\/4a. 
The junction conditions at 2: = give 

2(1 - 4ae)d,%, + (/? + 8afc)n47/.. = {^i^-'l {"If^- " ^) ^} (^.6) 

and 

Ah 

{l + (3k + Aak^)a^^ = ^S (3.7) 

where Sab = 2(5£mat/^^"* — /lab^mat IS the energy-momentum tensor for the brane matter. 

Switching to Fourier space, the bulk graviton equation ( |3.5| ) is solved, for spacelike 
momenta (p^ > 0), by 

7^,(p,z)oce2'=l^li^2(pe'=l^lA) . (3.8) 

We do not use the other solution as it diverges as 2; ^ 00. Substituting this into the 
junction conditions allows 7^,^, and hence the effective gravitational laws, to be determined 
on the brane. 

Define Q^u to be the Einstein tensor corresponding to the perturbation 7^4,^. At large 
distances {p <^ k) an asymptotic expansion of the solutions shows that, to leading order in 
w M^^S^_ii, with 

Ml = . 3.9 

k 

Therefore general relativity with effective four dimensional Planck mass Mpi is obtained. 

If at least one of a and P is non-zero, the linearised brane gravity at short distances 
{p 3> k) will instead be described by 

Q^u - 2(r?^^n4 - df,d^)^ ^ M-^S^,y (3.10) 

-204^ ^ m;^^s (3.11) 

where 

M|i = (/? + 8ak)M^ (3.12) 

2 _ 3(l + /3fc + 4afc^) , 

4(l-4aA;2) ' ^^'^^^ 

and (/9 oc 7. This is linearised Brans-Dicke gravity. M^'^ gives the coupling strength for 
the scalar degree of freedom. If Brans-Dicke theory is to be compatible with gravitational 
measurements of the solar system, we need MpjM^ < 10^^ ||5|. 

If a = /3 = the short distance gravity will be very non-standard, and this leads to 
strong constraints on the brane world parameters (1/A; < 0.1mm). In contrast, the short 
distance constraints on scalar-tensor gravity theories, and hence the above brane world, 
are quite weak (1/A; < 100 km if 4a/e^ is near to 1 ||ll|). 
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4. Einstein Branch 



A similar analysis to that in the previous section can be applied to the brane world solutions 
of section We use the perturbed metric ( |3.1| ), and take <I> = (l)Q—a\z\+if, with ip small. As 
before, and ip are determined by requiring that the brane remains at z = 0. The scalar 
field perturbations mix in with the trace of the metric perturbations, and this requires a 
slightly different choice of gauge. 

The 'Gauss-Bonnet' branches will be considered in the next section. We will begin with 
the 'Einstein' branch, which has k = —2(uj + l)a. We are interested in the parameter range 
< cj < — 1, and so in contrast to the Randall-Sundrum model, this solution will have 
negative warp factor. Not surprisingly the resulting brane gravity is significantly different 
to that in the previous section. 

It is useful to split 7^,^ into tensor and scalar parts 

1 4 /I dadA , 

li,v = lp.v^ -^lri^,y^ -\-r]^,y- ^-\X (4.1) 

with 

^ 8fcy - ai 

^ 4(2fc-a) ^ ' 

Requiring that the brane does not bend gives 

V' = -5.(x-<^) (4.3) 
a 

d^%v = . (4.5) 
Gauging away i?^, the remaining bulk equations (for z > 0) reduce to 

HO (dl - 2{2k - a)d, + e^'^'Dij = (4.6) 
Afioiuj + 4/3) (dl - 2{2k - a)d, + e^'^'aA x = (4.7) 



where fiQ = 1 — 4a(2(7 — k){a — k). Again there is the possibility of ghosts if the bulk 
curvature is too high. To avoid this we need fiQ > (and uj > —4/3). 

Taking i2mat to be independent of the brane junction conditions imply 

2^lo^,%. + (/? - 8a[2a - k]) 0,%, = {-^m- " ^ {v,u - ^) s| (4.8) 

2mo9.x + (/3 - 8a[2a - k]) □4X = " (2^ -^AP ^ ^^'^^ 

and 

(1 + [4a(2a - k) - 0][a - k]) a,{x - ^) = ^^J-^M^ ^ ^^"^"^ 
where = M^e-^"^"^. 
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h-a/k{-p/k) J 
- {2/io(\/<^2 + p2 + a) + (/? - 16aay | (4.14) 



For spacelike momenta the solution of the bulk equations (4^,4/7) which vanishes as 
z — > oo is 

e(2^-'^)^/2_./, (-pe^VA;) , (4.11) 
unless A: = (a; = —1), in which case it is 

exp a + 

Substituting in the bulk solution and setting z = 0, we find that the left hand side of 
the graviton equation ([4.81) is equal to 



for w > — 1 and 



for uj = —1. The left hand side of the scalar equation ( |4.9| ) will have a similar form. 

If the above term in brackets vanishes for some p > 0, then eqs. (^^) and (^^) will 
have non-trivial solutions when £mat = 0. In this case the vacuum will have tachyon modes, 
and so will be classically unstable. If the induced terms in the action (|2.3| ) are absent, or 
small {(3 < 16a [2 + w]), this is always the case. Hence the addition of the higher order 
gravity terms has destabilised the solution. 

If the quadratic curvature terms are absent, the above brane world has a negative 
tension brane, suggesting that the system is unstable. In fact it is stable, but only because 
of the Z2-symmetry. It is therefore not at all surprising that adding extra terms to the 
action (even if a is tiny) removes the stability. Similar effects were noted for the brane 
world models discussed in ref . ||l^ , and also for the two brane Randall-Sundrum model |^ . 
The problem could be avoided by taking /3 to be large. This will be discussed in section |6|. 

5. Gauss-Bonnet Branches 

We will now consider the new branches of solutions which appear when the quadratic 
curvature terms are included in the action. These have the bulk, background solution 



, 1 ujjuj + 2) 

We will only consider — 1 < u; < 0, since the field equations are degenerate when w = 0. As 
in the previous section 7^1, should be split into scalar and tensor parts. We use the same 
gauge as before (p!^4.5), but this time take 



1 ,4/1 a^a^ 



7/.^ = + -^IVf^u + ( -^Vfiu - ) X (5.2) 



where the field x is now defined by 



4(6k + c^a) 
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with 

Loa{3k + 2{2uj + l)cr) 



^ k{k + uja) 

Using the above expressions, the bulk equations reduce to 



(5.4) 



and 
where 



fi^ (a^ - 2{2k - a)d, + f^e^'^a^) = (5.5) 



di - 2{2k - a)d, + f'yu^ X = (5.6) 



2 k — {1 — uj)a 



' WO" + k 

" + 2{2u + l)a ^^'^^ 
H^ = 8a{k + 2[uj + l]a){k + uja) (5.9) 

fi^ = -16u;a{k + 2[u; + l]c7)(3A; + 2[2uj + l]cr) . (5.10) 

As in sections ^ and if the parameters /x^ and /Li;^ are not positive the solution will have 
bulk ghosts. 

In contrast to the Einstein gravity solutions of section ^ the tensor and scalar per- 
turbations (7^1/ and X respectively) have different behaviour in the bulk. For example, if 
/■y <C fx, then the tensor perturbations will be less localised on the brane than the scalar 
perturbations. In further contrast to the bulk field equations of the previous section, the 
parameters and may be negative. If this happens, the corresponding field equation 
will resemble that for a field in a spacetime with three timelike coordinates (or two if fi^^^ 
is also negative). This will give rise to instabilities on the quantum level and, for some 
parameter ranges, on the classical level too. Either way, we need to have not just /j.^ and 
fjL^ positive, but also and 

The brane junction conditions imply 

2fi^d,%u + + 8a[k - 2a])n,%, = j^^^ _ 1 (^^^^ _ s| (5.11) 

fix {2{k - a)d,x - fx^4x) + 2A; ([3 + 2uj]p + 8a[uj + 2] [3A: + 2uja]) □4^^ = -k^S (5.12) 

l^x^d^X - 3A;(/3 + I6a[k + uja])n4X + (4 + /3[3A: + 2uja])aa4ip = (5.13) 
with = M^e~'^'^° y/a. If the coefficient of the □47^;/ term in the boundary condi- 



tion ( 5.11| ) is negative, the system may be unstable. If A; < the situation will be similar 



to that in the previous section, and there will be a graviton tachyon mode. The same will 
be true in A; > and fx^ is not too small. Alternatively if /i^ is sufficiently small there 
will be no tachyon, but the effective four dimensional gravitational coupling will have the 
wrong sign. Similar considerations apply to the scalar modes. 

If we take the negative choice of sign for the solution (|5.lD , then the requirement that 
there are no ghosts restricts a to be less than cr* = [8a(2a;^ + 7uj + 6)]~^/^. In this range 
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8a[k — 2a] < and k is always negative. This branch is therefore unstable (for small /?). As 
with the solution in section Q the instability could be removed by including a large induced 
gravity term (i.e. taking /3 > —8a [A; — 2a]). We will consider this in section ^. 

For the rest of this section we will take the positive sign in eq. (^.1|). For this branch, 



k, fi-y, and are all positive, and > For to be positive we require that 
a < ac, where dc = [4a(u;^ + 2uj + 3)]~^^'^. To avoid tachyons (graviton and scalar) we also 
need f3 + 8a{k — 2a) > 0. If /? = this is true when o" < cr*. Since cr* < ac, the junction 
conditions provide a stronger constraint than the bulk effects. 
The bulk graviton equation (|5.5|) is solved by 

%,{p, z) oc e^^'-'^^K^^^/k {f^pe'^yk^ . (5.14) 

The bulk scalar equation has a similar solution, but with replaced by f^. The bulk 
graviton perturbations are qualitatively similar those in the conventional brane world (see 
section |3|). In fact we can get something similar to the Randall-Sundrum model (with a 
scalar field and induced gravity term) if we take a — > 0. In particular and will be 
both equal to one. 

The solution of the bulk field equation ( 5.141 ) implies that 



Ozlt,v = -{pf^)-^ ,,. l^lv (5.15) 

on the brane. A similar relation holds for dzX- By using asymptotic and series expansions of 
the above expression, we can determine the approximate small and large distance effective 
gravitational laws on the brane (just as we did in section ^). If A; > a (as is the case for 
the solutions we are considering), then when p is small dz^f_iu ~ ~f^p'^liMu/[2{k — a)], and 

when p is large d^^^u ~ -f-jPltiu- 

At large distances, when p <€.k/ f^, the junction conditions (^.11| - ^T3 ) imply that the 



four dimensional effective linearised gravity will be a tensor-scalar theory [see eqs. ( 3.1C| ) 
and (|ll|)] with 



Ml = [(3 + 2u;)(/3/V^) + 8V^(u; + 2)(3/c + 2uja)] . (5.17) 

In this case the effective brane scalar field, (f, is just (p. 

To avoid conflict with experimental measurements, we would like to have Mpi <C 
m. For id = —1 (so (7* = cjc) and (3 = this will be true when f-y « 0, which occurs 
when a ~ ac- We then have ~ and M-pi/M^ ~ Thus a hierarchy between 
the effective scalar and tensor couplings has arisen, despite the fact that the fundamental, 
five-dimensional scalar and gravitational couplings are of the same order. This is possible 
because the scalar field and the higher gravity terms have affected the different parts of 
the perturbation in different ways. This only occurs when the model includes both these 
features. If $ is turned off {a ^ 0), then — > 1 and both types of perturbation 'feel' 
the effects of the bulk in the same way. On the other hand if the higher gravity terms are 
dropped (a — > 0), the above solution does not exist in the first place. 
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A similar effect can be obtained for other values of uj, but unfortunately the region of 
parameter space with « (a « Uc) has 8a[k — 2a] < 0, which implies that tachyons 
will be present.. This can be fixed by including an appropriate induced gravity term with 
coefficient f3 = Pc = 8a(l + uj)ac- A hierarchy is then obtained, but at the price of a 
fine-tuning p. 

At very short distances, when p ^ k/ f^, the ^4 terms in the boundary conditions 
dominate the dz terms. The resulting effective gravitational theory is again scalar-tensor. 
The Planck mass is then 

M|i = + 8^{k - 2a)]M^ . (5.18) 

This time the effective scalar mode is a combination of 7 and (/?. This gives a different 
relation between the two mass scales. 



Ml = 3M|i 



^ ^Jlf \2Mf ^ ^ (Mpi/M*)2 + I6^{k + 2[1 + uj]a) 



(5.19) 

For the parameter ranges we are considering, this gives < 3Mpj, even when a is fine 
tuned to be near ac- The suppression of the scalar part of the gravity is lost at short 
distances. However the short distance constraints on Brans-Dicke gravity are weak, so this 
is not a problem if the length scale fx/^ is of geographical size or smaller. 

For k/ ^ p ^» k/ fx, tensor perturbations have large distance behaviour, while scalar 
perturbations have small distance behaviour. This third regime is another new feature of 
this brane world mode which is not present in the standard brane world. Mpi is again 
given by eq. (5.16) as above, but this time there is a different effective scalar, and 



M2A;((j-/c) / 3kY^ ^ ^' 

1+ 7T 2 + 



4:M^y/E{k + 2[l + u;]a)uja'^ \ a-k + ujgY 



(5.20) 



where 



MlAa-k) 

y = l + ^ ^ , '—— . (5.21) 



As with the short distance gravity, we find < 3Mpj. As before, this need not be a 
problem unless the length scale f-^/k is of astronomical size. 

The above expressions for Mpi and M^ at the three different scales are algebraically 
rather complicated, and not particularly revealing. It is instructive to take /3 = /3c, a = 
<7c — 1/0^) then consider e close to zero. In this case 



Ml, . "-^f^Mle (5.22) 
(-w) 



at large and medium distances, and 



M|,»?fci±fi^Af^ (5.23) 
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at short distances. At medium and short distances, it can be seen from eqs. ( 5.19 ) and 



( |5.20 ) that M? is approximately 3 times the corresponding Mpj (recah Mpi <C -M^ when 



f3 = f3c and a ^ CTc)- At large distances 

We see that Brans-Dicke gravity constraints could be satisfied (at large distance scales) by 
taking e < 10~^. On the other hand, we see from eqs. (|]22D and ( ^l23|) that there is at 



least a factor of 3 difference between Mpj on medium and small distance scales. This may 
lead to conflict with observations. Of course, any problems could be avoided by taking the 
distance scale /■y/k to be very small. In this case gravity on all scales at which accurate 
gravitational experiments have been performed would be given by eqs. ( 5.16| ) and ( [5.17 ). 



6. Scalar Gravity 

We saw in section ^ that if second order bulk curvature terms were added to a brane world 
solution with a negative warp factor, it could develop a tachyon mode. If a sufficiently 
large induced gravity term is added to the action, this instability will be removed. Having 
included such a term, the approximate gravitational law on the brane can be obtained (as 
with the other brane world solutions) by substituting the solutions to the linearised bulk 
equations (4.13, ^31 ) into the junction conditions (4.8,4.!i 



At relatively small distances {p ^> \k\), we obtain a scalar-tensor theory of grav- 
ity ( |3lO| , pl] ) with 

= [{|3/^/^) - 16{UJ + 2)V^C7] M2 (6.1) 



and 

Ml = {^ + 3u;) ,. • (6-2) 



(2 + a;)/Uo-fT(4 + 3w)M2j 



If w = —1 then M(f, = Mpi. For other values of lo it appears to be possible to obtain 
M(j) 3> Mpi by fine-tuning the denominator of eq. ( |6.2| ), but unfortunately M| is negative 
in this region of parameter space. Thus, for any choice of the above theory will not be 
suitable for describing gravity in our universe at astronomical distances. 

Looking at the behaviour of the perturbations for small momenta (large distances), 
we see that the junction conditions (4^,OD imply that j^u,X ~ while (for uj ^ — 1) 



99, 7 ~ S/p"^. To leading order in p the tensor modes are negligible and so at large distances 
(p <C |A;|), we have a scalar theory of gravity. A model with similar behaviour was discussed 
in ref. @]. 

By considering the motion of a test particle which is confined to the brane, we see that 
cPx^/dt^ ~ ~foO' ^° i^^^ static solutions) the gravitational potential is y = (1/2)700. 
For the above scalar theory of gravity there is an order 1/p'^ contribution to V from 7, and 
so a Newtonian potential is obtained. If on the other hand the particle is not confined to 
the brane, cPx^/dt'^ "^qq = — Fqq + kvi. Using the expression for Vi (|4.4|) , we find that 
now only 700 and x contribute to V, and so no Newtonian potential is obtained. Therefore 
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any Newtonian potential (at large distances) is purely a result of the brane's curvature, 
and nothing to do with the bulk behaviour of the perturbations. 

The large distance behaviour of the above model is also incompatible with astronomical 
observations. The addition of an induced gravity term to the solutions of section ^ may have 
solved the stability problem, but it still does not produce a theory which could represent 
our universe. If w = —1, we will instead have a non-standard, scalar-tensor, large distance 
gravitational law, with a non-Newtonian gravitational potential. We will leave the analysis 
of such a theory for future work. 

The 'Gauss-Bonnet' solutions in section ^ with the negative sign choice in eq. ( ^.1] ) 
suffered from a similar instability to the above 'Einstein' solution. This can also be cured 
by adding a large induced gravity term. In general, these stabilised solutions will also 
behave as scalar gravity theories with Newtonian gravitational potentials at large distances. 
At short distances they will be approximately Brans-Dicke, with Mpi and given by 
eqs. ( ^.18 ) and (|5.19 ). In contrast to the above 'Einstein' solution, it does appear to be 
possible to fine-tune /3 so that ^ Mpi without having < (although we will not 
discuss the algebraic details here). In this case a hierarchy is obtained by fine-tuning the 
positive contribution to the scalar coupling from the induced action (|2.3| ) , so that it cancels 
the negative contribution from the other boundary terms ( |2.2[ ). This contrasts with the 
hierarchy obtained in section ^, which arose from the properties of the bulk background 
solution. 

As with the other 'Gauss-Bonnet' solutions discussed in section ^, these solutions 
also have a medium distance gravitational law which is different from the small and large 
distance laws. Since this solution branch has < f^, the tensor modes will have short 
distance behaviour in this region, and the scalar modes will have large distance behaviour 
(the reverse of what happens for the other 'Gauss-Bonnet' solutions). Thus for k/f^^ 
p ^ k/ we will again have scalar-tensor gravity, but with different effective mass scales 
to the short distance theory. 



7. Conclusions 

In this paper we have investigated the effects of higher order curvature terms and a scalar 
field on linearised gravity in a brane world scenario. These two extensions of the con- 
ventional brane world give rise to many new features. The higher order gravity produces 
additional terms in the brane junction conditions. These resemble those produced by an 
induced gravity term, and allow gravity to be four dimensional at all length scales, weaken- 
ing constraints on the model. Another, less desirable, consequence is the possibility of new 
instabilities. If the curvature is too high, corrections to the bulk equations result in the 
appearance of ghosts. If the extra contribution to the junction conditions has the wrong 
sign, the model will have tachyons. 

Not surprisingly, the gravitational field equations for the quadratic theory have more 
solutions than the linear one. One of the two extra solution branches gives a better behaved 
brane world than the one obtained from the linear theory. In particular it is stable, and its 
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bulk space is warped in a similar way to the Randall-Sundrum model. The other solution 
branches are unstable when the higher order gravity effects are included. 

The introduction of the scalar field alters how perturbations 'feel' the warping of the 
bulk spacetime. The various coefficients in the linearised field equations are altered by 
presence of the bulk scalar field and the higher gravity terms. Significantly the changes 
to the scalar and tensor perturbation equations are not the same, and so with some fine- 
tuning of the background solution, it is possible to suppress the scalar modes relative to 
the tensor modes. This does not happen unless we include both the scalar field and the 
higher gravity terms. The effective four-dimensional theory on the brane is Brans-Dicke 
gravity which, thanks to the suppression of the scalar modes, could be compatible with 
observational constraints. However it should be noted that the constraints on Brans-Dicke 
gravity were obtained by going beyond a linearised approximation of the theory. Since 
we have only considered a linearised theory in this paper, the actual constraints may be 
different. 

The instabilities in the other solution branches can (in some cases) be removed by 
including sufficiently large induced gravity terms in the brane action. The resulting theory 
is then a standard four-dimensional gravity action, plus the brane world part of the model. 
Since Brans-Dicke theory is stable, it is not too surprising that the resulting combined 
theory will also be stable if the brane world part of it is small. The resulting brane gravity 
is approximately scalar-tensor at short distances. It appears that the scalar modes can be 
suppressed in some cases by tuning the size of the induced term to cancel the brane world 
contribution. A Newtonian potential is obtained at large distances, although (to leading 
order) the resulting theory has no tensor part. 

Having included quadratic curvature terms in the theory it is very natural to ask 
whether even higher order terms should be included as well. If the solutions in this pa- 
per had small curvature, we could safely ignore higher terms. Inevitably this is not the 
case, since if the curvature were small, the effects of the quadratic terms would also be 
small, and our brane world solutions would not be significantly different from the usual 
brane world scenarios with scalar fields. Work on brane worlds with higher order gravity 
actions [^] (without scalar fields) suggests that these higher curvature terms will give qual- 
itatively similar contributions to the linearised gravity as the quadratic ones. Hence the 
results obtained in this paper should remain qualitatively unchanged, although obviously 
the algebraic expressions for the different mass and length scales in the models will be 
different. 
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A. Appendix: Field Equations and Boundary Terms 



Variations of the various terms in the actions ( |2.1| ) and are given below. The results 
are derived for a D-dimensional manifold, M, bounded by a hypersurface S. The non- 
standard convention that the normal of S points into M. is used, so dPx^—g V^Fa = 
— Jj. d^~^x\/—h n°'Fa (for a brane world, this is slightly easier to visualise). Note that there 
is an extra factor of two in the brane world boundary terms, since the brane is treated as 
the boundaries of each of the two bulk half-spacetimes. In the expressions below ^ is an 
arbitrary function. In the rest of the paper it was equal to e~^*. 



5{ / d^x,/^iR-2 I d'^-'xV-h^K 



= [ d''xV^6g'''{^Gab-VaVb^ + 9abV^^} + 6nd^^R} 
Jm 

- [ d^'-'xV^dg'^'' {CiKab - Khab) - dn^hab} + 6^{2d^CK} 



(A.l) 



<5 / d^x^/^e(V$)2 
Jm 

= j^d^'x^ 6g''' jeVa W,«> - ^g,,{V^f^ + 5$ {a$e(V$)2 - 2V"eVa«' - 2CV2$} 
- j d^-^xyfH5^{2idn^] (A.2) 
Before dealing with the quadratic curvature terms, it is useful to define the tensors 



1 



Jab = ^[ 2KKacK% + K.dK^'Kab - 2KacK"'Kdb - K^Kab ) , (A.3) 



Pabcd — Rabcd + '^Rb[c3d]a " '^Ra[c3d]b + R9a[c9d]b 



SO £gb = P"''"^'^ Rabcd- The second order Lovelock Tensor M is then 



Pa^'^^Rbcde — -^gab^GB ■ 



1 



(A.4) 



(A.5) 



d^'x^^CcB- I d^^-^xV-hiaJ-'^G'^'Kab] 

M JT. 

[ d^x^ dg"^ {2iHab + 4P,,fee V^V^C} + 5^{d^iCGB} 
Jm 



d^-^xV-h6g^'i e ( a Jab - 2Jhab - ^PacbeK' 



r,o.b 



+ 2dni [^Gab + 2KeaK\ - 2KKab + hab[K^ " K.aK""]^ 
+ 4 {2K,^,D'^D,)^ - KDaDbi - K^bD^i - hab[K,,D'D% - KD^^]) | 
5^[Ad^S,{J -2G''''Kab)] (A.6) 
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+ 2K,(„L>'=$L>i,)$ - KDa^Dh^ - habK^'De^Dc^') 
+ dn^ [Df^a^Df,)^ - habD'^^Dc^] - ^On^ [Da^Db^ - KbiD^f] 

+ S^!^aK^-K''''Kab-R)dn^ 

+ (d^^D'^^D''^ - 2D"4L>^$ - 2CL>"£)^$) {Kab - Khab)^ (A.7) 



+ -(V$)2 [<7,,Ve?V^# - 2V(,eVb)<I>] 

(5$|2^[i?„5V"$V''$ + VaV6$V"V''$ - (V2$)2] + 4VaCV6$V"V''$ 

+ (5$^(V$)2 - 2V„^V"$)V2$ + v2,e(V$)^} 

- d^^dn^ (^(^"*)' + (^*)') - ^n^[(^n$)' + (A.8) 
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d^-^xV^5^ {4^dn<^> + } (A.9) 

References 

[1] K. Akama, An Early Proposal Of 'Brane World', Lect. Notes Phys. 176 (1982) 267 
[hep-th/0001113]; 

V. A. Rubakov and M. E. Shaposhnikov, Do We Live Inside A Domain Wall?, Phys. Lett. B 
125 (1983) 136; 

M. Visser, An Exotic Class Of Kaluza-Klein Models, Phys. Lett. B 159 (1985) 22 
[hep-th/9910093]; 

E. J. Squires, Dimensional Reduction Caused By A Cosmological Constant, Phys. Lett. B 
167 (1986) 286. 

[2] L. Randall and R. Sundrum, An alternative to compactification, Phys. Rev. Lett. 83 (1999) 
4690 [hep-th/9906064]. 

[3] J. Garriga and T. Tanaka, Gravity in the brane-world, Phys. Rev. Lett. 84 (2000) 2778 
[hep-th/9911055]. 

[4] O. De Wolfe, D. Z. Preedman, S. S. Gubser and A. Karch, Modeling the fifth dimension with 
scalars and gravity, Phys. Rev. D 62 (2000) 046008 [hcp-th/9909134]. 

C. Csaki, J. Erlich, T. J. Hollowood and Y. Shirman, Universal aspects of gravity localized on 
thick branes, Nucl. Phys. B 581 (2000) 309 [hep-th/0001033]. 

[5] T. Damour and K. Nordtvedt, Tensor-scalar cosmological models and their relaxation toward 
general relativity, Phys. Rev. D 48 (1993) 3436, (and references therein). 

[6] C. Lanczos, A Remarkable Property of the Riemann-Christoffel Tensor in Four Dimensions, 
Ann. Math. 39 (1938) 842. 

[7] D. Lovelock, The Einstein Tensor And Its Generalizations, J. Math. Phys. 12 (1971) 498. 

[8] J. E. Kim, B. Kyae and H. M. Lee, Effective Gauss-Bonnet interaction in Randall- Sundrum 

compactification, Phys. Rev. D 62 (2000) 045013 [hcp-ph/9912344]; 

Various modified solutions of the Randall- Sundrum model with the Gauss-Bonnet interaction, 
Nucl. Phys. B 582 (2000) 296 (Erratum-ibid. B 591 (2000) 587) [hcp-th/0004005]. 
J. E. Kim and H. M. Lee, Gravity in the Einstein-Gauss-Bonnet theory with the 
Randall- Sundrum background, Nucl. Phys. B 602 (2001) 346 (Erratum-ibid. B 619 (2001) 
763) [hcp-th/0010093]. 

[9] I. P. Neupane, Consistency of higher derivative gravity in the brane background, JHEP 0009 
(2000) 040 [hep-th/0008190]. 

Y. M. Cho, I. P. Neupane and P. S. Wesson, No ghost state of Gauss-Bonnet interaction in 
warped background, Nucl. Phys. B 621 (2002) 388 [hcp-th/0104227]. 

I. P. Neupane, Completely localized gravity with higher curvature terms. Class. Quant. Grav. 
19 (2002) 5507 [hep-th/0106100]. 

Y. M. Cho and I. P. Neupane, Warped brane-world compactification with Gauss-Bonnet term. 
Int. J. Mod. Phys. A 18 (2003) 2703 [hcp-th/0112227]. 

[10] K. A. Meissncr and M. Olcchowski, Domain walls without cosmological constant in higher 
order gravity, Phys. Rev. Lett. 86 (2001) 3708 [hep-th/0009122]; 



-17- 



Brane localization of gravity in higher derivative theory, Phys. Rev. D 65 (2002) 064017 
[hep-th/0106203]. 

[11] N. Deruelle and M. Sasaki, Newton's law on an Einstein 'Gauss-Bonnet' brane, Prog. Theor. 
Phys. 110 (2003) 441 [gr-qc/0306032]. 

[12] C. Charmousis, S. C. Davis and J. F. Dufaux, Scalar brane backgrounds in higher order 
curvature gravity, JHEP 0312 (2003) 029 [hep-th/0309083]. 

[13] S. C. Davis, Generalised Israel junction conditions for a Gauss-Bonnet brane world, Phys. 
Rev. D 67 (2003) 024030 [hcp-th/0208205]. 

E. Gravanis and S. Willison, Israel conditions for the Gauss-Bonnet theory and the 
Friedmann equation on the brane universe, Phys. Lett. B 562 (2003) 118 [hep-th/0209076]. 

[14] G. W. Gibbons and S. W. Hawking, Action Integrals And Partition Functions In Quantum 
Gravity, Phys. Rev. D 15 (1977) 2752. 

[15] R. C. Myers, Higher Derivative Gravity, Surface Terms And String Theory, Phys. Rev. D 36 
(1987) 392. 

[16] N. E. Mavromatos and J. Rizos, String inspired higher- curvature terms and the 
Randall-Sundrum scenario, Phys. Rev. D 62 (2000) 124004 [hep-th/0008074]; 
Exact solutions and the cosmological constant problem in dilatonic domain wall 
higher-curvature string gravity. Int. J. Mod. Phys. A 18 (2003) 57 [hep-th/0205299]. 

[17] A. Jakobek, K. A. Meissner and M. Olechowski, New brane solutions in higher order gravity, 
Nucl. Phys. B 645 (2002) 217 [hep-th/0206254]. 

[18] P. Binetruy, C. Charmousis, S. C. Davis and .J. F. Dufaux, Avoidance of naked singularities 
in dilatonic brane world scenarios with a Gauss-Bonnet term, Phys. Lett. B 544 (2002) 183 
[hep-th/0206089]. 

[19] I. Y. Aref'eva, M. G. Ivanov, W. Muck, K. S. Viswanathan and I. V. Volovich, Consistent 
linearized gravity in brane backgrounds, Nucl. Phys. B 590 (2000) 273 [hep-th/0004114]. 

[20] H. CoUins and B. Holdom, Linearized gravity about a brane, Phys. Rev. D 62 (2000) 124008 
[hep-th/0006158]. 

[21] C. Charmousis and J. F. Dufaux, Gauss-Bonnet gravity renders negative tension branewolds 
unstable, hep-th/0311267. 



-18- 



